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ABSTRACT 

In this paper, we describe an analytical method for treating uniformly rotating ho- 
mogeneous rings without a central body in Newtonian gravity. We employ series ex- 
pansions about the thin ring limit and use the fact that in this limit the cross-section 
of the ring tends to a circle. The coefficients can in principle be determined up to 
an arbitrary order. Results are presented here to the 20th order and compared with 
numerical results. 
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1 INTRODUCTION 

The problem of the self-gravitating r ing captured the in - 
terest of such reno wned scientis t s as iKowalewskvl l|l885l ). 
iPoincare! (|l885l l and lDvsonl (|l892|jl893h . Each of them tack- 
led the problem of an axially symmetric, homogeneous ring 
in equilibrium by expanding it about the thin ring limit. In 
particular, Dyson provided a solution to fourth order in the 
parameter a — a/b, where a provides a measure for the ra- 
dius of the cross-section of the ring and b the distance of the 
cross-section's centre of mass from the a xis of rotation . 

First numerical results were given by I Won g| (|l974l ). who 
was not able to clarify t h e tran sition to spheroidal bod- 
ies. | Eriguchi fe Sugimotol (Il98ll ) and lEriguchi fc Hachisul 
(1985) developed improved methods with which they 
were able to study this transition and the connection 
to the Maclaurin spheroids. Returning to the p roblem 
significantly later, lAnsorg. Kleinwachter fc Meinell (|2003(J ) 
achieved near-machine accuracy, which allowed them to 
study bifurcation sequences in detail and correct erroneous 
results. 

In this paper we provide a scheme to extend Dyson's 
work up to arbitrary order in a. With the help of com- 
puter algebr4E we were able to compute the solution ex- 
plicitly up to the 20th order and see a marked improvement 
over the fourth order. Although constant mass density is 
of considerable importance to the method presented here, 
an extension of it to other equations of state is possible 
|Petroff fc HoratscheklliuOg ). 



* E-mail: S . HoratschekStpi .uni-jena.de (SH); 
D.Petroff@tpi.uni-jena.de (DP) 

1 We made use of Maple™ . Maple is a trademark of Waterloo 
Maple Inc. 
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Figure 1. A sketch providing the meaning of the coordinates 
(r, x) in relation to the cylindrical coordinates (g, z). 



2 THE APPROXIMATION SCHEME 
2.1 The Coordinates 

To describe axially symmetric rings, we introduce the polar- 
like coordinates (r, \, p), which are related to the cylindrical 
coordinates (g, z, tp) by 



g = b — rcosx, z = rsm\, (p = tp, 



(1) 



see also Fig. [TJ The surface of the ring will be described by 
a function r s (x)- In our coordinates, the Laplace operator 
applied to a function / = f(r, x) reads 



_ j. d 2 f 1 df 1 d' 2 f 
V / = — - H H — 

dr 2 r dr r 2 d\ 2 



— (b — r cos x) I cos x 



df sin x df 



dr 



(2) 



We choose the constant b in such a manner, that the centre 
of mass of the ring's cross-section coincide with r = 0, thus 
implying 



2tt r a (x) 



/ir cos^dr dx = 0, 



(3) 



where a is the mass density. 
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2.2 Basic Equations 

For solving the problem of a self-gravitating fluid in equilib- 
rium, we have to fulfil Laplace's equation 



V 2 U = 

outside the fluid and Poisson's equation 
V 2 U = 



(4) 



(5) 



inside it, where U is the gravitational potential. Addition- 
ally, we have to satisfy Euler's equation 



dv „ TT _ 

M "d7 = " M ~ p ' 



(6) 



where v is the velocity of a fluid element and p the pressure. 
We consider uniform rotation about the axis g — with the 
angular velocity fl and thus have the velocity field 



v = ft x x 



leading to 



V \ U + 



dp' 



Integration gives 

U + 



f dp' 1 
J M(P') 2 



(7) 



(8) 



(9) 



where Vq is the constant of integration. At the surface of the 
ring, where the pressure vanishes, we have 



= V . 



(10) 



We are dealing here with homogeneous rings, i.e. the mass 
density is a constant 



fi = constant. 



(11) 



This means that the integral in Q and is simply p/n 
and ((3]) reduces to 



^7T 

/ 



r s (x) cosxdx = 0. 



(12) 



2.3 The Idea of the Approximation Scheme 

The thin ring limit is approached when the ratio of the inner 
radius gi to the outer one g a tends to 1. In this limit, the 
cross-section of the ring becomes a circle. This is the start- 
ing point of the approximation. As an interesting aside, if 
one considers a ring surrounding a central body, for exam- 
ple a point mass, the cross-section of the ring can deviate 
significantly fr om a circle even if the ra dius ratio is close to 
1, see Fig. 8 in lAnsorg fc Petrofj (|2005l ). 

We describe the surface of the ring by the Fourier series 

r s (x)^a(l + J2J2l3 lk co S (kx)a l + o(a q )) , (13) 

\ i=l fe=l / 

where 

a:=°-. (14) 



There are no sine terms because of refiectional symmetry 
with respect to the equatori al plane, which is k nown to hold 
for fluids in equilibrium fsee lLichtensteinll 19331 ). To the lead- 
ing order, the cross-section is indeed a circle of radius a. We 
make a similar ansatz for the square of the angular velocity 



n 2 = ttGa* (^2 + °( aq+2 )j 



and Vo 



■kG^o 2 ( Via 1 + o(a q ) 



(15) 



(16) 



Because we consider constant mass density, a given sur- 
face function r a (x) determines the distribution of the mass 
completely, i.e. one can compute the potential everywhere. 
For technical reasons, we first calculate the potential on the 
axis of symmetry only. After this, we determine the potential 
outside the ring and in particular along its surface, where 
(TTOl) must hold. 

Since the potential outside the ring contains logarithmic 
terms in r it will come as no surprise that there are In a terms 
in the coefficients of our series. Like Dyson, we introduce 



a 



(17) 



Let us suppose that we already know the solution up to 
the (q— l)st order, and want to find the qth order, in partic- 
ular f3 q k (fc = 1, 2, . . . , q) and fi 9 +i. To find these q + 1 un- 
knowns, we have to fulfil equation (|12p and also require that 
the q Fourier coefficients in front of cos(fcx) (fc = 1> 2, • • • , q) 
in (|10|) vanish. After finding these q + 1 unknowns, we now 
know Vq-i, but cannot yet determine the coefficient v q . 



2.4 The Potential Outside the Ring 

First off, the Poisson integral is used to find the value of the 
potential along the axis of rotation. We label the coordinates 
for a point on the axis (R,Xr) — ( r >x)i from which 



b = Rcos(xr) 
follows. The axis potential is 

2tt r s ( x ') 

U^ ia {R) = -2nG 

o o 

2vrr s ( x ') 



(18) 



/j,(b — r cos x') r 



s/R?T: 



2Rr cos %l> 



dr dx 



= — 2-kG^i f I (b — r cosx) ("j?) Pi{cos^)drdx 



-2tt G^a^a a ) — i R2l l 



(19) 



where ip := — xr- Please note that the expansion in terms 
of powers of 1/R indicated in the last line is not trivial, 
since there is an independence hidden in the terms with ip. 
Because of refiectional symmetry, there are only terms with 
odd powers in 1/R. We expand Ai with respect to a 



(20) 
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Using (Taxis, we can then find the potential anywhere in the 
vacuum region. To do so, we first introduce a set of axially 
symmetric solutions to Laplace's equation that vanish at 
infinity. We define 



h(g,z) 



7T 

/ 



dip 



z 2 — 2bgcos ip 



(21) 



which is nothing other than a multiple of the potential of a 
circular line of mass with radius b, centred around the axis. 
In (r, x, (p)-coordinates it reads 



2K 



h(r, X ) 



4fr 2 — 4br cos x \ 
4b' 2 — 4br cos x + r ' 2 ) 

Ab 2 — Abr cos \ + r2 



(22) 



where K denotes the complete elliptic integral of the first 
kind, 



K(k) 



_' 



(10 



v/l - k 2 ! 



(23) 



Because the difference of two of such solutions with different 
6's also satisfies Laplace's equation, it is clear that 

l-l 

(24) 



1 d 



where 







db db ^ dr 



d sinx d 



(25) 



r dx 

is also a solution of Laplace's equation. Next we note that 
along the axis we have 

ty(21 - 1)!! 



(26) 



(21 - l)R 2l ~ l ' 

It then follows that the potential in the vacuum region is 
/?+! 

U out (r,x) = -2-KGfia 2 ^ a 21 ' 1 ^ 1 Aih{r, X ) + o(a q ) 



(27) 

since this expression satisfies Laplace's equation, vanishes at 
infinity and has the correct value along the axis. For calcu- 
lating the potential at the body's surface we expand Ii (r, x) 
for r < b. For example we get 



h(r,x) 



.M»,~ 

r 



ln(8fc/r) — 1] cosx t 



+ 



•0 



(28) 



After evaluating these equations for Ii (r, x) at the surface 
r — r B ( X ), we use (|7T) l to find the coefficients 4>ih in the 
expansion 

(7 s (x) = -2-KGiia 2 J2 4>ik casikxW + o{a q ) J . (29) 



2.5 The Potential Inside the Ring and the 
Pressure 

To treat the interior of the ring, it is convenient to introduce 
the dimensionless coordinate 

y ■= -• (30) 



For calculating the potential inside the ring, we first expand 
it, 

q i 

Uin(r,x) = -ivGiia 2 ^2^2U ik {y) cos(k X )cr l ■ (31) 

«=0 k=0 

Inserting this into Poisson's equation ((5} and collecting the 
coefficients in front of cos(kx)o^ ^ gives ordinary differential 
equations (ODEs) in y for the Uik{y)- The solutions to these 
equations are determined uniquely by requiring that Uik{y) 
be finite at the centre y — and that the potential be contin- 
uous at the surface, i.e. Ui a (r B (x),x) must agree with U B (x)- 
After solving the ODEs for Uik(y), we know the inner 
potential, and with <(9j we are able to calculate the pressure 



r, x) = nGfi 2 a 2 ^2 ^Pik{y) cos(fcx)cr l + °{° q *) 



(32) 



3 SOLUTION UP TO THE FIRST ORDER 

To this order, we shall show the calculation in detail. The 
unknown quantities are /3n, f2o, Qi, Q2 and vo- The surface 
function reads 



r s (x) = o[l + cos x + o(<t)] 
and equation (|12p becomes 

2-k 



(33) 



/ 



r s (x) 3 cos X dx = a 3 [3n finer + o(a)]. (34) 



This gives 

fin = 0, (35) 

cf. Table lA"2l which means that we have the surface function 

r.(x) = o[l + o(<7)]. (36) 

Next we have to evaluate equation (|10fl • We immediately see 
that 



fio = = 
must hold, and furthermore we have 
1 



(37) 



- 2 2 

2« e 



ivGfia 



— +[ — - Q 2 cos x ) cr + o(a) 



(38) 



at the surface. To calculate the potential at the ring's surface 
to the same order in a, we first compute the potential alon 
the axis. Here we have to calculate the potential of a torui 



87rG/ia 3 
3R 



(39) 



The expansion in terms of powers of 1/R gives 

Ai = l + o(<r), (40) 

2 denotes the complete elliptic integral of the second kind, 
E(k) := f ( f \f\ - fe 2 sin 2 9d6». 
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and 



thus 



A 2 = --a + o(a) 



Qio = 1, an = and 021 = 



8' 



(41) 



(42) 



cf. Table IA4I Now we are able to calculate the potential 
at the body's surface via (I27|l . (|28|l and the corresponding 
equation for I 2 . Using equation (|36p it is possible to expand 
the surface potential in a. We get 



U S ( X ) = -2TvGfia 

which implies 

4>oo = A + 2, 



and 



(43) 



(44) 



Plugging JTSJl, (J38j) and (g3j) into equation (|T0j) and col- 
lecting the coefficients in a l cos(fcx), we find the following 
equations: 

i k equation 



2A + 4+±fi 2 

1 ±fi 3 = vi 

1 1 A + § - n 2 = 

Solving these equations gives 



: Co 



cf. TableHU and 



n 2 



VO 



A+ ■ 



}A + 



(45) 



(46) 



The equation with O3 and vi cannot be further evaluated 
until the next order q — 2. 

For the mass and the angular momentum, we get 



M = 2ttVV~ 1 +o(ct -1 )] 



and 



J = y/ir 5 Gfj. 3 (4A + 3) a 5 [a' 2 + o(a~ 2 )]. 



(47) 



(48) 



To leading order, Qi/g = 1 — 2<r holds, and we can conclude 
that 



lim 

Pi/Po — 



47rbVo 



5GM 



In 1 



go 



In 16, 



(49) 



see a lso equation (11) in iFischer. Horatschek fc Ansord 
l|2005l ). A comparison of this result for homogeneous rings 
with the analogue for poly tropic rings can be found in 
iPetroff fc Horatschek! <|2008h . 

To calculate the inner potential, we have to find a so- 
lution to Poisson's equation. The ansatz (I3ip leads to the 
ODE 

d 2 U 00 (y) , ldC/ 00 (y) 



+ 



dy 2 ' y Ay 
to leading order, which has the solution 



+ 4 = 



(50) 



At the centre, the potential has to be regular, thus Ci — 0. 
The resulting potential at the surface r = r s (x) is 



U s (x) = irG(Ma 2 ll-C 2 + o(l) 



(52) 



For the potential to be continuous, C2 = 2A + 5 must hold, 
see (|43fl . To the zeroth order the potential is not a function 
of the angle \- 

Furthermore we have the equations 



d 2 U w (y) | l dUi (y) _ Q 



dy 2 



V dy 



and 



d 2 Un(y) , ldUu(y) U u (y) 



dy 2 



+ - 



y dy 



y 



+ 2y = 0. 



(53) 



(54) 



Note that the 2y term in the second equation results from 
Uoo (y) , which is already known. The solutions of these equa- 
tions that are regular at the centre and have the correct 
values at the surface are 



and 



cf. Table 



and 



Uio(y) = 



Uu(y) = (X + l)y-^, 
For the pressure, (J9J leads to 

Poo(y) = l-y 2 , 
Pw(y) = 



(55) 



(56) 



(57) 
(58) 



(59) 



After finding the inner potential and pressure, we can cal- 
culate the potential energy 



W := 2 J UdV = -^G^a 5 
the rotational energy 



T := 



g 2 dV = 7r 3 G/i 2 a' J 



2A+ H ) a-' + oia- 1 ) 



X+^a-'+oia- 1 ) 



and the integral over the pressure 

P := j pdV = 7r 3 GViV (a' 1 +o(a- 1 )) (60) 



to first order. We see that the virial theorem (|62p is fulfilled 
up to this order. 



4 DISCUSSION 

With this approximation method, we are able to calculate 
e.g. the shape, angular velocity and pressure of the ring up 
to arbitrary order in a. We have done so up to the 20th 
order. 

To test our solutions, we ensured that the transition 
condition 



Uoo(y) 



-y +Ci my + C 2 . 



(51) 



VIU S = Vt/outl 



(61) 
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Table 1. For a given radius ratio gi/g = 0.9, physical quantities to different orders in q and numerically 
determined values are compared to the values for q = 20: M 2 o = 4.6299179884304816293 X 10"" 2 , S°2| = 
3.2474683264953211610 xl0~ 2 , J 20 = 7.5456215256289320669 x 10" 3 , P 20 = 1.7862946528142761708 x 10~ 4 , 
T 20 = 6.7988816964653749490 x 10" 4 , W 20 = -1.8956647351373578410 x 10" 3 . 





<j 


Mg/M 20 


- 1 




- 1 


Jq/ J20 ~ 


1 


Pq/P20 ~ 


1 


fg/f 20 ~ 


1 


W q /W 20 


- 1 


1 


0.053 


1 x 10" 


-2 


1 x 10 


-2 


2 X 10" 


-2 


3 x 10" 


-2 


2 x 10" 


-2 


2 x 10" 


-2 


2 


0.052 


2 x 10" 


-4 


6 x 10 


-4 


4 x 10" 


-4 


6 x 10" 


-3 


-1 x 10" 


-3 


6 x 10" 


-4 


3 


0.052 


2 x 10" 


-4 


1 x 10" 


-4 


2 x 10" 


-4 


4 x 10" 


-4 


3 x 10 


-4 


3 x 10 


-4 


■1 


0.052 


3 x 10" 


-6 


-1 x 10 


-5 


-3 x 10" 


-6 


5 x 10" 


-5 


-5 x 10" 


-5 


-2 x 10" 


-5 


10 


0.052 


6 x 10" 


11 


-9 x 10" 


11 


4 x 10- 


12 


5 x 10- 


10 


-6 x 10~ 


10 


-3 x 10- 


10 


19 


0.052 


4 x 10" 


17 


4 x 10- 


17 


6 x 10- 


17 


1 x 10- 


16 


9 x 10- 


17 


9 x 10- 


17 


num 




-4 x 10" 


14 


1 x 10" 


16 


-4 x 10 _ 


14 


-7 x 10" 


14 


-4 x 10- 


14 


-3 x 10" 


14 



Table 2. For a given radius ratio Q\/g = 0.5, physical quantities to different orders in q are compared 
to the numerically determined values M num = 0.7201292, Q 2 um = 0.5467604, J num = 0.3247949, P num = 
0.04874713, T num = 0.1200820, W num = -0.3864053. 





(7 


Mg/M num 


- 1 


i£ num 


- 1 


Jq/ ^num 


- 1 


Pq 1 Pnum 


- 1 


Pq /Tnum 


- 1 


Wg/Wnum 


- 1 


1 


0.33 


2.8 x 10" 


1 


2.3 x 10- 


1 


4.2 x 10 


-1 


8.6 x 10 


-1 


4.6 x 10 


-1 


6.1 x 10" 


1 


2 


0.30 


5.1 x 10- 


-2 


5.9 x 10" 


-2 


6.5 x 10 


-2 


2.5 x 10 


-1 


2.6 x 10 


-2 


1.1 x 10- 


-1 


3 


0.30 


5.1 x 10- 


-2 


4.2 x 10- 


-2 


5.8 x 10 


-2 


1.3 x 10 


-1 


8.1 x 10 


-2 


1.0 x 10- 


-1 


4 


0.30 


1.7 x 10" 


-2 


1.3 x 10- 


-2 


1.9 x 10 


-2 


5.4 x 10 


-2 


8.4 x 10 


-3 


2.6 x 10- 


-2 


10 


0.29 


1.2 x 10- 


-3 


1.0 x 10- 


-3 


1.3 x 10 


-3 


3.2 x 10 


-3 


1.0 x 10 


-3 


1.9 x 10" 


-3 


20 


0.29 


2.6 x 10" 


-5 


2.5 x 10" 


-5 


3.0 x 10 


-5 


7.2 x 10 


-5 


2.7 x 10 


-5 


4.4 x 10" 


-5 



is fulfilled up to the appropriate order in a. Furthermore we 
tested that the virial theorem 



W + 2T + 3P = 



(62) 



is fulfilled for each order in a. 

An important question is how good this method is. In 
Tables [l] [2] and [3] one can see how the dimensionless quan- 
tities 



M 
JI 



1 ' 

Wo 



P 
P 



f 

' T 



1 

- YL 

~ w ' 

with 



J 
J 



1 



C-l/2^3/2^,5' 



(63) 



G/i 2 £»o 

increasing order for differ- 



miprove m accuracy 
ent radius ratios. Especially for thin rings, we get very 
accurate results. In fact, for rings with radius ratios 
Qi/ Qo ~ 0.85 we achieve a precision which is compara- 
ble with that given by the numerical m ethod described in 
lAnsorg. Kleinwachter fc Meinell (l2003al ). For larger radius 
ratios, the accuracy is thus better. As a co-product, our 
work provides an independent test of the accuracy of the 
numerical method (better than 10" 13 cf. Table [TJ. 

The shape of the ring in meridional cross-section for 
various radius ratios can be found in Fig. [2] The curves to 
order q — 20 can barely be distinguished from the numerical 
ones for Qi/g > 0.3. As one approaches the transition to 
spheroidal topologies (qi/Qo — » 0), the true curve becomes 
pointy at the inner edge and is no longer well represented by 
our Fourier series. Nevertheless, the shape of the ring is quite 
well approximated even for Qi/ g a = 0.1, as seen in Fig. [3] 
The surface function r s (x)> which is a constant to leading 
order, clearly approaches the numerical one with increasing 
q. The pressure in the equatorial plane can also be seen 
to approach the numerically determined one for Qi/g = 
0.3 in Fig. [4] It is interesting to note that the centre of 
mass does not coincide with the point of maximum pressure. 
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Figure 2. Meridional cross-sections of rings to the order q = 
20 for different radius ratios Qi/g - The g- and z-axis are scaled 
identically in such a manner that g Q has the same value for all the 
rings. The dot in each ring marks the centre of mass of the cross- 
section (g = b, z = i.e. r = 0) and the dashed line shows the 
numerical result and is indistinguishable from the q = 20 curve 
for Qi/g > 0.3. 



In Fig. [5] one can get an impression of the accuracy of the 
approximation over the whole range of radius ratios and for 
various values of q. Despite the claim found in IWonel (|l974l ) 
that Dyson's perturb ative metho d diverges for a > 1/3 (see 
also the comments in lDvsonlll892j) . these results indicate the 
opposite. 

The (dimensionless) coefficients 0^ and Q,i are only 
functions of the dimensionless quantity a = a/6, but not 
of any length like q for example. The reason for this is that 
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Table 3. For a given radius ratio Qi/Qo = 0.2, physical quantities to different orders in q are compared 
to the numerically determined values A? num = 0.9424, fi£ um = 0.9844, J num = 0.4545, P n um = 0.07865, 
f num = 0.2255, Wnum = -0.6869. 



q 


a 


Mg/M num 


- 1 


"?/"num 


- 1 


Jq 1 ^num " 


- 1 


Pq 1 ^num 1 


Tq /Tnum 


- 1 


W q /Wnum 


- 1 


i 


0.67 


1.0 




7.5 x 10 


-i 


1.6 




5.1 


1.6 




2.8 




2 


0.54 


3.2 x 10" 


l 


3.5 x 10 


-l 


3.9 x 10" 


-l 


1.6 


3.2 x 10 


-l 


7.6 x 10" 


l 


3 


0.54 


3.2 x 10" 


l 


2.7 x 10 


-l 


3.5 x 10" 


-l 


1.1 


5.4 x 10 


-l 


7.2 x 10" 


l 


4 


0.51 


1.9 x 10" 


l 


1.7 x 10 


-l 


2.1 x 10" 


-l 


6.3 x lO" 1 


2.1 x 10 


-l 


3.6 x 10" 


l 


10 


0.47 


6.4 x 10" 


2 


6.9 x 10 


-2 


7.1 x 10" 


-2 


1.9 x 10" 1 


8.7 x 10 


-2 


1.2 x 10" 


l 


20 


0.46 


2.2 x 10" 


2 


2.7 x 10 


-2 


2.3 x 10" 


-2 


6.0 x 10~ 2 


3.2 x 10 


-2 


4.2 x 10" 


2 




r s (x) 
b 



0.9 
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X 

Figure 3. The meridional cross-section and the (dimensionless) 
surface function r B (x)/b of the ring with radius ratio Qi/t}o = 
0.1 for different orders q compared to the numerical result. The 
surfaces are scaled such that g a (and therefore also £;) has the 
same value to all orders. 
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Figure 4. The pressure in the equatorial plane for a ring with 
radius ratio Si/ Qo = 0.3 for different orders q compared to the 
numerical result. It is interesting to note that the centre of mass 
of the cross-section does not coincide with the point of maximum 
pressure. To the order q = 18 we get b := (6— Bi)/(@o — Bi) = 0.503 
and p := (g p ,max — 8\)/ (<?o ~ Qi) = 0.480, which differ by less than 
1% from the numerical values. 



equations @, (JE), © and are scale invariant in the 
following sense: If U(x), fi(x), p(x) and v(x) are solutions 
of these equations then a~ 2 U(ax), fj,(ax), a~ 2 p(ax) and 
a^ 1 v(ax) are too, where a is an arbitrary scaling factor. 
Other equations of state fi — fi(p) will not satisfy this scal- 
ing invariance in general. To parametrize a specific ring, in- 
cluding its dimension, we need two parameters (for example 
a and 6) as expected. 

Please note that one has to be careful in interpreting 
the results for the thin ring limit. For example, one might 
think that the squared angular velocity vanishes like a 2 In a. 
This is true for the dimensionless quantity Q 2 /Gfj,, but need 
not be true for the squared angular velocity itself. If we fix 
the 'size' b and the mass M of the ring in that limit, then 
the cross-section shrinks to a point (a — ab). With (|47l) we 
can conclude that fj, oc a~ 2 and therefore Q 2 cx ln<r, which 
means that Q, 2 and hence the velocity of a fluid element go 
to infinity. 

Relativistic rings, including the thin ring limit , 
were studied in | Ansorg. Kleinwachter fc Meinell ( 2003bl ). 
lAnsorg et all (|2004T i and lFischer et all (120051 ). From the per- 
spective of General Relativity, the Newtonian theory con- 
stitutes a good approximation when certain conditions are 
fulfilled. For one thing, typical velocities must be small com- 



pared to the speed of light and for another G\U\ <? must 
hold. We just saw, however, that for rings of finite extent 
and mass, the velocities grow unboundedly in the thin ring 
limit. The same holds for U s as well, see (152 p . This means 
that the Newtonian theory of gravity is not appropriate to 
describe this subtle limit itself, since one cannot expect it 
to be a good approximation to General Relativity. It is re- 
markable that the approximation about the point a — is 
nevertheless so successful. 
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Table Al. Coefficients Qi up to the order q = 9 (i 
Hi 



2,3,. 



2 
3 
1 
5 
6 
7 
8 
9 

10 



_1 \ _ 19 
8 96 



25 X 3 i 365 \2 _ 2345 \ _ 8989 
128 ~ r 1536 18432 73728 



25 >4 , 40235 \3 _ 134255 x 2 _ 36493505 x _ 34831813 
64 A ~r 98304 A 294912 A 56623104 A 226492416 



134925 \5 , 2797535 x 4 _ 50072105 \3 _ 
131072 1572864 169869312 

1021727845 \2 _ 104581877693 \ _ 49377918425 
509607936 97844723712 391378894848 



APPENDIX A: FURTHER COEFFICIENTS 

In Tables [ATUA4I coefficients for fli, flik, Uik(y) and au are 
given. 



Table A2. Coefficients f3ik up to the order q = 9. The bold-faced type indicates that these terms are incorrect in lDvsonl 1 1892T ) 
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1325 \3 I 
8192 " r 
109915 \2 i 
442368 "r 

863633 , 
18874368 

105833029 
4076863488 
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Table A3. Coefficients Uik(y) up to the order q = 7. 
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Table A4. Coefficients an up to the order q = 8. 
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